arXiv: 1509.0421 Ivl [cs.IT] 14 Sep 2015 


Wireless Throughput and Energy Efficiency with 
Random Arrivals and Statistical Queueing Constraints 

Mustafa Ozmen and M. Cenk Gursoy 


Abstract —Throughput and energy efficiency in fading channels 
are studied in the presence of randomly arriving data and statistical 
queueing constraints. In particular, Markovian arrival models in¬ 
cluding discrete-time Markov, Markov fluid, and Markov-modulated 
Poisson sources are considered. Employing the effective bandwidth 
of time-varying sources and effective capacity of time-varying wire¬ 
less transmissions, maximum average arrival rates in the presence of 
statistical queueing constraints are characterized. For the two-state 
(ON/OFF) source models, throughput is determined in closed-form 
as a function of the source statistics, channel characteristics, and 
quality of service (QoS) constraints. Throughput is further studied 
in certain asymptotic regimes. Furthermore, energy efficiency is 
analyzed by determining the minimum energy per bit and wideband 
slope in the low signal-to-nolse ratio (SNR) regime. Overall, the 
impact of source characteristics, QoS requirements, and channel 
fading correlations on the throughput and energy efficiency of 
wireless systems is Identified. 

Index Terms —effective bandwidth, effective capacity, energy 
efficiency, fading channels, Markovian source models, maximum 
average arrival rates, minimum energy per bit, QoS provisioning, 
wideband slope, wireless throughput. 

I. Introduction 

A. Motivation 

Mobile data traffic has experienced unprecedented growth 
recently and is predicted to grow even further over the coming 
years. For instance, it is projected that global mobile data traffic, 
which already grew 81 percent in 2013, will increase 11-fold 
between 2013 and 2018, reaching 15.9 exabytes per month by 
2018 Cl. As one of the main driving forces behind this growth, 
wireless transmission of multimedia content has significantly 
increased in volume and is expected to be the dominant traffic in 
data communications. Indeed, mobile video traffic was already 
53 percent of the traffic by the end of 2013 and is predicted 
to increase 14-fold between 2013 and 2018, accounting for over 
two-thirds of the world’s mobile data traffic by 2018 Cl. 

This exponential growth in the flow of mobile data and multi- 
media content has significant implications on wireless networks. 
For one, wireless multimedia traffic requires certain quality-of- 
service (QoS) guarantees. For instance, in voice over IP (VoIP), 
multimedia streaming, interactive video, and online gaming appli¬ 
cations, constraints on delay, packet loss, or buffer overflow prob¬ 
abilities need to be imposed so that acceptable performance and 
quality levels can be met for the end-users. Another consequence 
is heterogeneity in network traffic. Wireless networks now carry 


heterogeneous traffic in diverse environments, and successful 
design of networks, efficient use of resources, and effective QoS 
provisioning for multimedia communications critically depend on 
the appropriate choice of source traffic models. For instance, 
while voice traffic can be accurately modeled as an ON/OFF 
process with fixed-rate data arrivals in the ON state, data traffic 
can be bursty and video traffic, which exhibits correlations, 
can be modeled statistically using autoregressive, Markovian, or 
Markov-modulated processes m. 

Finally, it is important to note that this increased traffic 
together with the given QoS requirements need to be supported 
by wireless systems equipped with only limited bandwidth and 
power resources. Especially, due to limited energy available for 
mobile units and rising energy costs and environmental concerns, 
energy efficiency in wireless communications is a key concern 
(see e.g., |[3l and H). Therefore, it is crucial to identify the 
fundamental performance limits (e.g., in terms of maximum 
achievable throughput and minimum energy per bit) in order to 
determine how to most effectively utilize the scarce resources. 
With this motivation, in this paper we investigate the throughput 
and energy efficiency in fading channels when data arrivals are 
in general random, and QoS constraints in the form of limitations 
on the asymptotic buffer overflow probabilities are imposed. 

B. Literature Overview 

Satisfying QoS requirements is critical for most communica¬ 
tion networks, and how to satisfy QoS constraints for various 
source traffic models has been one of the key considerations 
in the networking literature. In particular, besides conventional 
queueing theory, network calculus has been introduced by Cruz 
in early 1990s as a theory to address the delay and other 
deterministic service guarantees in networks by dealing with 
queueing systems a-Ei. Subsequently, Chang in lH developed 
the stochastic version of the network calculus. More specifically, 
the theory of effective bandwidth of a time-varying source has 
been formulated to identify the minimum amount of transmission 
rate that is needed to satisfy the statistical QoS requirements 
(see also 13 - mi This theory is based on the logarithmic 
moment generating function of the arrival process and is re¬ 
lated to the large deviation principle. Moreover, statistical QoS 
constraints are imposed as limitations on buffer/delay violation 
probabilities. Effective bandwidths of various source models 
have been investigated extensively in the literature. Eor instance, 
Elwalid and Mitra studied the effective bandwidth of Markovian 
traffic sources (including Markov-modulated fluid and Markov- 
modulated Poisson sources) in IH under constraints on the 
buffer overflow probability. It is shown that effective bandwidth 
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is given by the maximum eigenvalue of a matrix derived from 
source parameters and service requirements. In iTsll . effective 
bandwidth formulations were provided for multi-class Markov 
fluids as well as memoryless (Poisson) and discrete-time Markov 
sources. In El, the authors studied the effective bandwidths of 
general stationary sources and derived a hrst order approximation 
of the effective bandwidth in terms of the mean arrival rate and 
index of dispersion. 

In wireless communications, the instantaneous channel capac¬ 
ity varies randomly depending on the channel conditions. Hence, 
in addition to the source characteristics, transmission rates for 
reliable communication are also time-varying. In such cases, 
randomly time-varying servers can be considered in the queueing 
system model. Indeed, motivated by the wireless channel, Stol- 
yar in El, Venkataramanan and Lin in El, and Sadiq and 
de Veciana in El employed tools from the theory of large 
deviations and investigated scheduling rules (e.g., MaxWeight, 
Exponential, and Radial Sum-Rate Monotonic scheduling) while 
controlling the large deviations of queues. Following another 
method, the time-varying channel capacity can be incorporated 
into the theory of effective bandwidth by regarding the channel 
service process as a time-varying source with negative rate and 
using the source multiplexing rule (' ifTTl Example 9.2.2]). Using 
a similar approach, as a dual concept to effective bandwidth, Wu 
and Negi defined in 1201 the effective capacity, which describes 
the maximum constant arrival rate that a given time-varying 
service process can support while satisfying the statistical QoS 
requirements. Indeed, work in l20l revitalized the consideration 
of statistical queueing constraints in the context of wireless com¬ 
munications, and the effective capacity of wireless transmissions 
has been investigated intensively in various settings (see e.g., 
im-Ga). For instance, Tang and Zhang in l22l considered 
the effective capacity when both the receiver and transmitter 
know the instantaneous channel gains, and derived the optimal 
power policy that maximizes the system throughput under QoS 
constraints. Liu et al. in 12^ considered fixed-rate transmis¬ 
sion schemes and analyzed the effective capacity and related 
resource requirements for Markov wireless channel models and 
Markov fluid sources. In l26l and El, effective capacity of 
cognitive radio channels was studied. In l28l . multi-antenna 
communication in the presence of queueing limitations was 
investigated. Soret et al. in l30l addressed correlated Rayleigh 
fading channels and studied the effective capacity under different 
adaptive rate policies. In this study, performance in the presence 
of probabilistic delay constraints and variable rate sources was 
also analyzed by considering a Gaussian autoregressive source 
model. Energy efficiency in the presence of QoS limitations was 
addressed in El - ina. 


C. Contributions 

We note that the studies on the effective capacity of wireless 
channels have primarily concentrated on constant arrival rates in 


the analysis of the throughput and energy efficiencjQ. Departing 
from this approach, we in this paper explicitly take into account 
the randomness and burstiness of the source traffic. In particular, 
we address Markovian source models including discrete-time 
Markov, Markov fluid, and Markov modulated Poisson sources, 
and conduct a performance analysis. More specihcally, our con¬ 
tributions can be listed as follows: 

• A framework with which source randomness can be incor¬ 
porated in the throughput analysis of wireless transmissions 
is provided. 

• For two-state (ON/OFF) source models, closed-form expres¬ 
sions are obtained for the maximum average arrival rate in 
terms of the source statistics, effective capacity of wireless 
transmissions, and the QoS exponent 6, which quantihes 
how strict the QoS constraints are. 

• Throughput is characterized in the low-0 and high-SNR 
regimes. 

• An energy efficiency analysis is conducted and minimum 
energy per bit and wideband slope expressions are deter¬ 
mined for both constant and random arrival models. 

• Via both analytical and numerical results, the impact of 
source randomness, fading correlations, and queueing con¬ 
straints on the wireless throughput and energy efficiency is 
identified. 

The remainder of this paper is organized as follows. In Section 
III-AI we describe the channel model. Sections III-BI and III-CI 
contain the preliminaries regarding the statistical queueing con¬ 
straints, effective bandwidth, and effective capacity. In Section 
imi we provide our characterizations of the throughput with 
Markovian source models by analyzing the maximum average 
arrival rates. We conduct an energy efficiency analysis in Section 
|IV]for both constant and Markovian arrivals. Finally, concluding 
remarks are given in Section [V] Proofs are relegated to the 
Appendix. 

IT System Model 

A. Channel Model 

As depicted in Fig. [T] we consider a point-to-point link with 
a single transmitter and single receiver. In this system, the 
data generated by the source is initially stored in a buffer at 
the transmitter before it is transmitted over a wireless channel. 
We consider a flat-fading channel between the transmitter and 
receiver, and assume a block-fading model with a block duration 
of m symbols. Hence, fading varies independently from one 
block to another. On the other hand, we further assume that 
in each block duration of m symbols, fading can be arbitrarily 
correlated. The channel input-output relation within each block 

'To the best of our knowledge, the two exceptions to this are references 
(23] and (30) as also described in Section ITb] However, these studies have 
different modeling assumptions for the sources a nd/o r wireless transmissions 
from what we have in this paper. For instance, in (23), while ON/OEF Markov 
fluid anivals are considered, wireless transmissions occur at fixed rates and 
wireless link is als o m odeled as a continuous-time Markov chain with ON 
and OFF states. In a Gaussian autoregressive source is considered in the 
analysis. Additionally, our analysis, maximum average aiTival rate expressions, 
and throughput characterizations in the low-0, high- and low-SNR regimes are 
novel contributions with respect to these prior studies. 
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Fig. 1. System Model. 


can be expressed as 

Hi = h^Xi +ni for i = 1, 2,..., m (1) 


where Xi and yi are the channel input and output, respectively. 
The average energy of the input is E, i.e., 

E{|a:,p}=f. (2) 


Ui denotes the zero-mean, circularly-symmetric, complex Gaus¬ 
sian noise with variance E{|nip} = Nq. Hence, the signal-to- 
noise ratio is 


SNR = 


E{|:rp} 




(3) 


E{|nP} 

Above in O, hi denotes the fading coefficient. Fading coeffi¬ 
cients are assumed to be identically distributed, and the fading 
distribution can be arbitrary with finite variance. 

While the ensuing analysis is applicable to a general class 
of fading distributions, we use a Gauss-Markov fading model 
in the numerical results and assume that the Gaussian fading 
coefficients in each block of m symbols follow the correlation 
pattern hi = p/ii-i +Wi where Wi is an independent, zero-mean 
Gaussian random variable with variance E{|uiip} = (1 — 
p G [0,1], and is the common variance of the fading 
coefficients {hi}. Note that when p = 1, we have full correlation, 
whereas p = 0 models the case of independent fading. 


B. Queueing Constraints 


We assume that the data to be transmitted is generated from 
random sources and is first stored in a buffer before transmission. 
Statistical constraints are imposed on the queue length. In par¬ 
ticular, we assume that the buffer violation/overflow probability 
satisfies 


lim = -9 


(4) 


(j->oo q 


where Q denotes the stationary queue length, and 9 is the decay 
rate of the tail distribution of the queue length. The above limiting 
formula implies that for large q, we have 


MQ >q}- e-®«. (5) 


Indeed, a closer approximation is ll20l 

Pr{(5 > g} « (6) 

where = Pr{(5 > 0} is the probability of non-empty buffeid 
From (|6]l, we notice that, for a sufficiently large threshold, the 


^Probability of non-empty buffer can be approximated from the ratio of average 
arrival rate to average service rate [33] 


buffer overflow probability should decay exponentially with rate 
controlled by the QoS exponent 9. Note that as 9 increases, 
stricter queueing or QoS constraints are imposed, while looser 
queueing constraints are implied by smaller values of 9. Con¬ 
versely, for a given buffer threshold q and overflow probability 
limit e = PrjQ > q}, the desired value of 9 can be determined 
as 

6»=iloge-. (7) 

q e 

In the given setting, the delay violation probability is also 
characterized to decay exponentially and is approximated by ED 

Pr{£) > 4 « (8) 

where D is the queueing delay in the buffer at steady state, d is 
the delay threshold, and a* (9) is the effective bandwidth of the 
arrival process, described below. 

Next, we introduce the notions of effective bandwidth and 
effective capacity which we subsequently employ to formulate 
the wireless throughput in fading channels in the presence of 
random arrivals and statistical queueing constraints. 

1) Effective Bandwidth: Effective bandwidth characterizes 

the minimum constant transmission (or service) rate required 
to support the given random data arrival process while the 
buffer overflow probability is limited or more explicitly the 
statistical queueing constraint described by (|4|l is satisfied. Let 
{a{k),k = 1,2,...} be a sequence of nonnegative random 
variables, describing the random arrival rates. Let also the time- 
accumulated arrival process be denoted by A{t) — CL{k). 

Then, the effective bandwidth is given by the asymptotic loga¬ 
rithmic moment generating function of A{t) E), i.e., 

a*(4= lim ilogE|e®^(‘)|. (9) 

t-J-oo Ot I J 

In Section III-CI we describe the effective bandwidth of differ¬ 
ent source arrival models in detail. 

2 ) Effective Capacity: Effective capacity, as a dual concept to 

effective bandwidth, identifies the maximum constant arrival rate 
that can be supported by a given time-varying service process 
while satisfying (|4]i. Let {i^[fc], fc = 1, 2,...} denote the discrete¬ 
time stationary and ergodic stochastic service process and S[t] = 
J2k=i '^he time-accumulated service process. Then, the 

effective capacity is given by ll20ll 

^^(SNR,^) = - lim ^log,E(e-®^[‘l|. (10) 

t->oo tfl t J 

Note that we have assumed that the fading coefficients {hi} 
change independently from one block of m symbols to another. 
Under this assumption, effective capacity simplifies to 

C£;(SNR,4 = -ilog,E{e-®"} (11) 

where i/ is the instantaneous service (or equivalently transmis¬ 
sion) rate in one block. If the channel input sequence {xi} is 
an independent and identically distributed (i.i.d.) sequence of 
Gaussian random variables with zero mean and variance £, then 
the service rate can be written as 

m 

log2(l + SNRZj) (12) 

i=l 
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where we have defined Zi = \hi\’^. Hence, the effective capacity 
in the units of bits/block is 

Cis(SNR,6») =. (13) 


and pi 2 = l—pii - For these transition probabilities, we can easily 
see that the probability of the ON state in the steady state is 


F^on 


1 -Pii 

2-pii -P 22 ' 


( 20 ) 


Remark 1: In the special case of independent channel coef¬ 
ficients in each block and Rayleigh fading, we can express the 
effective capacity in closed-form as 


m . 


Cis(SNR,6»)= -Trlog, 


SNR i°Se 2 eSNR r 1 - 


logg2 ’ SNR 


TTl Tfi 

= m log 2 (SNR) - ^ - J log, r( 1 - 


(14) 
1 


log ,2 ’ SNR 

(15) 


where r(s, tu) = dr is the upper incomplete gamma 

function. 


Therefore, the average arrival rate is 


favg — 


AFon = A 


1 -Pii 

2 -Pii -P 22 


( 21 ) 


which is equal to the average departure rate when the queue is 
in steady state ||39]| . 

2) Markov Fluid Sources: In this subsection, we address 
Markov fluid sources where the source arrival process is modeled 
as a continuous-time Markov chain. Assume that G is the 
irreducible transition rate matrix of the Markov chain, Ai is the 
arrival rate in the state, and A = diag {Ai, A 2 ,..., A„}. Then, 
the effective bandwidth of this source is given by m, M 


C. Effective Bandwidths of Different Source Models 

1) Discrete-Time Markov Sources: In this subsection, we 
consider discrete-time Markov source models. Assume that the 
transition probability matrix of the n-state irreducible and ape¬ 
riodic Markov source process is denoted by J, and Ai is the 
arrival rate in state i. Moreover, A = diag {Ai, A 2 ,..., A„} is 
the diagonal matrix of arrival rates. Then, the effective bandwidth 
of this discrete Markov source is given by im 

a(^) = ^loge [sp (e^'^J)] (16) 

where sp(-) is the spectral radius of the input matrix. Note that 
the stationary distribution tt can be found from the solution of 

ttI = 1, 

ttJ = 7r (17) 


A + ic) (22) 

where /r(-) denotes the maximum real eigenvalue of the input 
matrix. We also note that the stationary distribution tt of the 
continuous-time Markov chain can be found by solving 

7rl = 1, 

ttG = 0 (23) 

where tt = [tti, 7r2,..., 7r„], 0 = [0,...,0]^ and 1 = 

In order to derive closed-form expressions in our analysis, we 
again consider two states (ON/OFF). When there is no arrival, the 
state is OFF. When the state is ON, the arrival rate is A bits/block. 
The transition rate matrix for a two-state Markov fluid is in the 
form of 



where tt = [tti, 7r2,..., 7r„] and 1 = [1,..., 1]"^. 

In order to unveil the key relationships and tradeoffs, we 
consider a particularly simple two-state model. We assume that 
data arrival is either in the ON or OFF state in each block duration 
of m symbols. When the state is ON, A bits arrive (i.e., the arrival 
rate is A bits/block), while there are no arrivals in the OFF state. 
For this two-state model, the state transition probability matrix 
is given as 


J = 


Pii 

P 21 


P 12 

P 22 


(18) 


Given the above transition matrix J, the effective bandwidth for 
this ON-OFF Markov model can be derived as m 


G = 




(24) 


where a is the transition rate from OFF state to ON state whereas 
P is the transition rate from ON state to OFF state. Using (l22l l. 
we can express the effective bandwidth as 


a*{e) = 4 - (a + /3) + {OX - {a + P)Y + AadX 

lu L 


(25) 


The probability of ON state, 7r2, is required to define the average 
rate. Inserting the generator matrix G in (l24l l into (l2Tt . we obtain 
the ON state probability as 


1,/^ X-I li (pil+P22e^^ + \J (pii+P22e^®)2-4(pii+p22-l)e-'^ \ 

(p, A) — g loggi 2 I 

(19) 


where pn denotes the probability of staying in the OFF state 
from one block to another. Similarly, P 22 denotes the probability 
of staying in the ON state. The probabilities of transitioning from 
one state to a different one are therefore denoted by p 2 i = 1—P 22 


7J‘2 = FoN = 


a 

a + P' 


(26) 


Therefore, the average arrival rate of the two-state Markov fluid 
process is 


ravg — 


APon = A 


a 

a-\- P 


(27) 
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3) Markov Modulated Poisson Sources: In this subsection, we 
assume that the data arrival to the buffer is a Poisson process 
whose intensity is controlled by a continuous-time Markov chain. 
For instance, the intensity of the Poisson arrival process is Ai in 
the state of the Markov chain. Therefore, the source arrival 
is modeled as a Markov-modulated Poisson process (MMPP). 
Assuming that the G is the irreducible transition rate matrix of 
the Markov chain and A = diag {Ai, A 2 ,..., An} is the diagonal 
matrix of the intensities of the Poisson arrivals in different states, 
the effective bandwidth is given by m, m 

a*(0) = ip((e®-l)A + G). (28) 

As in previous sections, we consider a two-state (ON/OFF) 
model in which there are no arrivals in the OFF state (i.e., the 
intensity is 0) and the intensity of the Poisson arrival process is A 
bits/block in the ON state. Assuming the same generator matrix 
G as in (l24l i. we can express the effective bandwidth as 

[(e«-l)A-(a + ^)] 

+ -l)X-{a + ^)]^ + Aa (e« - 1) A. (29) 

Note that the average arrival rate in bits/block is again given by 
ravg = APon = A ——. (30) 

a + p 

We further note that if the transition rate (3 = 0, then we have 
PoN = 1. In this case, MMPP model specializes to a pure Poisson 
source with intensity A bits/block, and the effective bandwidth 
of this source is given by 

(31) 

III. Throughput with Markovian Source Models 

In this section, we formulate the throughput of wireless 
fading channels when the data arrivals are random and sta¬ 
tistical queueing constraints are imposed. More specifically, 
we consider Markovian arrival models introduced in Section 
III-Cl namely discrete-time Markov sources, Markov fluids and 
Markov-modulated Poisson arrivals. The states in these Markov 
processes are differentiated by the corresponding arrival rates in 
these states, e.g., the arrival rate in the state is Xi. If the 
stationary distribution of the Markov process is denoted by tt, 
the average arrival rate in an n-state Markov source model simply 
becomes 

n 

Tavg — ^ ^ '^iXi (32) 

i=l 

which is equal to the average departure rate when the queue is 
in steady state lf39l . 

We seek to determine the throughput by identifying the max¬ 
imum average arrival rate that can be supported by the fading 
channel described in Section III-AI while satisfying the statistical 
QoS limitations given in the form in (|4]i. As shown in ll^ Theo¬ 
rem 2.1], dH is satisfied, i.e., buffer violation probability decays 
exponentially fast with rate controlled by the QoS exponent 9, 


if the effective bandwidth of the arrival process is equal to the 
effective capacity of the service process, i.e., 

a*{9) =Ce{s^t^,9). (33) 

Hence by solving (l33l l. we can determine the maximum average 
arrival rate r*^,g(SNR, 6). By specifying the effective bandwidth of 
different source models and incorporating the effective capacity 
of time-varying wireless transmissions in Cl, the maximum 
average arrival rate can be determined for general n-state Marko¬ 
vian source models. Indeed, several n-state source models are 
addressed in Section |IV] However, in our analysis in this section, 
to illustrate the impact of the arrival and system parameters in a 
lucid setting, we concentrate on the two-state (ON-OFF) arrival 
models and provide closed-form expressions for the maximum 
average arrival rates in terms of the source parameters and the 
effective capacity of the wireless transmissions. We also identify 
the characteristics of the throughput in the low-6* and high-SNR 
regimes. We note that the analysis throughout this section is 
applicable to any arbitrary fading correlation within each fading 
block, with the exception of high-SNR characterizations which 
are obtained under the assumption of i.i.d. fading. 


A. Discrete-Time Markov Sources 


In this section, we consider two-state (ON/OFF) discrete 
Markov sources described in Section III-Cll and initially char¬ 
acterize the maximum average arrival rate that can be 
supported by the fading channel while satisfying the statistical 
QoS limitations given in the form in (|4]i. 

Theorem 1: For the two-state (ON/OFF) discrete Markov 
source, the maximum average arrival rate (in bits/block) as 
a function of the QoS exponent 9, effective capacity of the 
fading channel ((/^(SNR, 0), and the state transition probabilities 
is expressed as 


r:,g(SNR, 9) 



g2ecE(SNR,e) _ p^^gecE(SNR,e) 

1 - Pit - P22 + 


(34) 


Proof: See Appendix lAl 

Note that r* above is formulated in terms of the effective 

avg 

capacity, Ce, of wireless transmissions. In Fig. |2l we plot the 
the maximum average arrival rate as a function of the effective 
capacity for different source characteristics when 9 = 1. It is 
easy to verify that when Pqn = 2 -p^i-p 22 = 1 or equivalently 
P 22 — 1, (1^ simplifies to r}^,g(SNR,0) = C£:(snr, 0). Hence, 
when the source is always ON and therefore the arrivals are at 
a constant rate, maximum average arrival rate is equal to the 
effective capacity, as also observed in Fig. |2] On the other hand, 
we notice in this figure that as Pqn diminishes and the source 
becomes more bursty, throughput diminishes as well and smaller 
average arrival rates are supported for given effective capacity. 

As also indicated in the above discussion and seen in (l3^ . 
r}^g(SNR, 0) is in general a function of the state transition 
probabilities of the Markov arrival process in the presence of 
buffer constraints. On the other hand, as shown in the following 
result, this dependence disappears if no buffer constraints are 
imposed, i.e., when 0 = 0. 
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Fig. 2. Maximum average arrival rate rJvg vs. effective capacity C£;(SNR) for 
different source statistics. No fading correlation, i.e., p = 0. 0 = 1. 


distributed with unit mean (i.e., Rayleigh fading is experienced). 
Then, we have 


Pqn 


log2 SNR + 0(1) if > 


-r:^g(SNR,6») = 


Slog: 

PoN l0g2 SNR + 0(1) 
log2 SNR + 0(1) 


log2 
if 0 < 6» < 

if 6» = 0 


log2 e 


(38) 


as SNR —>■ 00 . 

Proof: See Appendix O 

Note that the high-SNR slope is defined as ifJTl 


5. 


lim 

SNR—^oo 


—r* „(sNR, 9) 

m avg V 7 / 

log2 SNR 


(39) 


Theorem |3] shows that the high-SNR slope of the maximum 
arrival rate for the two-state discrete Markov source that can be 
supported in the i.i.d Rayleigh fading channel is 


Pqn 

e log 2 e 
^ON 
1 


log2 

if 0 < 6» < 
if 6> = 0 


e 

1 


loga e 


(40) 


Theorem 2: As the statistical queueing constraints are relaxed 
by letting the QoS exponent 0 approach zero, the maximum 
average arrival rate converges to 

m 

lim r*^,g(SNR, 9) — E{log 2 (l + SNRz^)} bits/block. (35) 

i=l 

Moreover, the first derivative of with respect to 0 at 0 = 0 
is 


9r* (SNR, 9) 


39 


0=0 


= - - y]] COv{log2 (1 + SNRZj ), log2 ( 1 + SNRZj ) } 


* j'=l 


1 

2 


y^E{log2(l + SNRZ,)} j (36) 


Vi=l 


where we define p as 


9 = 


(1 -P 22 )(P 11 +P 22 ) 

(1 - Pii)(2 - pii - P 22 ) 


(37) 


It is interesting to observe from Theorem |2] that when no buffer 
constraints are imposed i.e., when 9 — 0, the high-SNR slope 
is Soo — 1. again independent of the source statistics. On the 
other hand, when 0 > 0, Soo becomes proportional to the ON 
probability and is now less than one unless the arrival rate is 
constant. Furthermore, for 9 values greater than -, Soo starts 
decreasing with increasing 9. Hence, the result in Theorem [3 
quantifies the performance degradation experienced at high SNR 
levels due to source randomness and statistical buffer constraints. 

Let us further simplify the source model and set pn = 1 — s 
and P 22 = s. The source is now described by the single parameter 
s. Notice that with this choice we have Pqn = s and hence s 
becomes a measure of the burstiness of the source. The smaller 
the s, the less frequently the data arrives and the more bursty the 
source becomes. At the other extreme, if s = 1, source is ON all 
the time and we have constant arrival rate. Furthermore, with the 
above choice of pu and P 22 , the expression for the maximum 
average arrival rate simplifies to 


Proof: See Appendix iBl 

We see from that if no statistical buffer constraints are 
imposed i.e., if 0 = 0, then the maximum average arrival rate 
is equal to the ergodic capacity of the block-fading channel, and 
therefore is independent of the statistical characteristics of the 
discrete Markov arrival model. Moreover, the dependence of the 
maximum arrival rate in this regime on the channel statistics is 
only through the marginal distributions of the fading coefficients. 
Hence, channel correlation in each fading block does not play 
any role. However, this radically changes when 0 > 0. For 
instance, we notice from (l36l l that even with a small increase 
in 0, r*^,g starts varying with the source and channel statistics, as 
exemplified by the dependence of the first derivative on p and 
the covariance function. 

Having discussed the low-0 regime above, we next provide 
a characterization of r*^,g (SNR, 0) at high SNR values for i.i.d. 
Rayleigh fading. 

Theorem 3: Assume that the channel fading coefficients are 
i.i.d. in each block and fading power z = is exponentially 


Cg(SNR, 0) 


■loge 


,6»CE(SNR,e) _ n - 


(1-s) 


(41) 


which can readily be seen to be a diminishing function as 
s decreases. Therefore, source burstiness generally hurts the 
throughput if we keep all other variables fixed. 

We can further observe this in Fig. |3 where we plot the 
maximum average arrival rate (or equivalently the throughput) as 
a function of SNR for different values of s and the QoS exponent 
0. Numerical analysis verifies that as the source becomes more 
bursty with lower values of s, throughput diminishes. Conversely, 
throughput is maximized when s = 1 i.e., when we have constant 
arrival rates. It is also interesting to notice from drill that the 

1 • I I • I • ■ 1 \ « PvelSNR.e) 

arrival rate in the ON state, which is given by A = -, 

increases as s diminishes. Hence, smaller s implies that data 
arrives less frequently but with bursts of increased rates. We also 
observe in Fig. [3 that the throughput reduction due to burstiness 
is more severe at high SNRs. This is indeed a consequence of 
the fact that high-SNR slope gets smaller as Pqn = s decreases. 
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Fig. 3. Maximum average arrival rate vs. signal-to-noise ratio SNR for 
different values of 6 and source statistics. No fading coirelation, i.e., p = 0. 



Fig. 4. Required SNR vs. ON probability, Pqn. for a given fixed average arrival 
rate. No fading correlation, i.e., p = 0. 


as discussed above. Finally, we see in Fig. [3 that performance 
degradation is experienced as 0 increases and hence stricter buffer 
constraints are imposed. 

In Fig. m we plot the SNR levels required to support a given 
average arrival rate as a function of the ON-state probability for 
different values of the QoS exponent 0. We observe that as Pqn 
decreases and hence the source becomes more bursty, required 
SNR level increases in general. Interestingly, a sharper increase is 
experienced under stricter buffer constraints (e.g., when 6 — 0.5 
rather than 0 = 0.1), indicating higher power/energy costs in 
these cases. 

The low-0 regime is investigated in Fig. |5] where we plot 
the maximum average arrival rate vs. QoS exponent 6 for 
different Pqn values. We set SNR = 1. We notice that all three 
curves converge to the same throughput value r*^g(0) as 0 —?> 0, 
confirming the result in (l35T l. Hence, source characteristics do not 
affect the throughput if no queuing constraints are imposed. As 
0 increases, throughput diminishes and the reduction in r* is 



Fig. 5. Maximum average arrival rate vs. QoS exponent 9 for different 
values of PoN- No fading correlation, i.e., p = 0. 



Fig. 6. Maximum average arrival rate r*,g vs. QoS exponent 9 for different 
fading correlations. Pqn = 0.5. 


more severe for more bursty sources (e.g., when Pqn = 0.4). We 
notice that, as predicted by (l3^ . this is already reflected by the 
different slopes of in the vicinity of 0 = 0. Hence, overall 
the system for more bursty sources becomes more cautious and 
supports smaller average arrival rates in order to avoid buffer 
overflows. 

In Fig. |6] we again plot the throughput as a function 0 but for 
different values of p, which quantifies the correlation between 
fading coefficients in each fading block. We fix SNR = 1 and set 
^ON = 0.5. Similar to burstiness, fading correlation does not have 
any effect on the throughput when 0 = 0. When 0 > 0, higher 
correlation (i.e., larger p) results in lower supported throughput 
under the same QoS constraints. 

Finally, we have conducted simulations to further verify the 
theoretical analysis and results. In particular, in the simulations, 
for fixed QoS exponent 0, SNR, and state transition probabilities 
Pii and P 22 of the ON/OFF discrete Markov source, we initially 
determine the maximum average arrival rate from (l34l l and the 
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Fig. 7. Buffer overflow probability PrjQ > g} vs. buffer threshold q for 
different values of d. pn = P 22 = 0.8, SNR = 0 dB 


Fig. 8. Delay violation probability Pr{D > d} vs. delay threshold d for 
different values of 9. pn = p 22 = 0.8, SNR = 0 dB 


corresponding maximum arrival rate in the ON state. Then, 
using the given statistical characterizations and the maximum 
arrival rate, we generate random Markov arrivals and assume 
that the arriving data is initially stored in the buffer before 
being transmitted. Transmission rates are simulated by generating 
realizations of i.i.d. Gaussian fading coefficients. Throughout this 
process, we track the queue evolution and the buffer state (i.e., 
the queue length) as the Markov arrivals occur (and hence more 
data gets stored) and transmissions at varying rates according to 
the generated fading coefficients are performed, clearing some 
data off the buffer. In Figs. [T] and [HI we plot the simulated 
buffer overflow probability Pr{Q > q} and delay violation 
probability Pr{£) > d}, respectively, as functions of the cor¬ 
responding thresholds, following 10^ runs of the simulation. We 
notice that while the theoretical analysis makes use of results 
from the theory of large deviations and is generally applicable 
for large thresholds, the simulation results are interestingly in 
excellent agreement with the theoretical predictions even at small 
values of the thresholds. For instance, we note from (jh]) that 
loggPr{Q > (?} ~ loge? — dq and hence is expected to decay 
linearly in q with slope 9. We indeed observe this linear decay in 
Fig. [7] (where the overflow probabilities are plotted in logarithmic 
scale) for even small to moderate values of q. Moreover, the 
slopes of the simulated curves, denoted by dsim, are very close 
to the originally selected value of 9. Similar conclusions apply 
to Fig. HI as well. In this figure, delay violation probabilities are 
determined by keeping track of the delay experienced by the 
data stored in the buffer until transmission. We again notice that 
the logarithm of the delay violation probability decays linearly 
with threshold d (or equivalently the delay violation probability 
diminishes exponentially with d). Note that the slope of the 
linear decay is predicted from ([H]) to be 9a* (9) where a* (9) is 
the effective bandwidth of the source. Again, the slope of the 
simulated curves are almost the same as this theoretical slope 
value, as indicated in the legend on the figure. 


B. Markov Fluid Sources 


In this section, we consider Markov fluid sources. In the fol¬ 
lowing, we go through similar steps as in the previous subsection 
and initially determine the maximum average arrival rates of 
ON/OFF Markov fluid sources that can be supported by the 
wireless channel as a function of the source transition rates and 
the effective capacity of wireless transmissions. Subsequently, we 
give characterizations of the maximum average arrival rates in the 
low-0 and high-SNR regimes. 

Theorem 4: For the two-state (ON/OFF) Markov fluid source, 
the maximum average arrival rate is given as 


r;,g(SNR, 9) = PoN 


^(^^(sNR, 9) -\- a + fi 
9 Ce{s^^, 9) + a 


(^^(SNR,^). 


(42) 


Proof: See Appendix iDl 

Note that maximum average arrival rate generally depends on 
the transition rate matrix of the Markov fluid source. At the 
same time, similar to the discrete case, when there are no QoS 
constraints, source characteristics do not have any impact on 
the throughput. However, this changes drastically when 0 > 0 
even if 0 is vanishingly small. These properties are demonstrated 
analytically in the result below. 

Theorem 5: As the statistical queueing constraints are relaxed 
by letting the QoS exponent 0 approach zero, we have 


m 

lim (SNR, 0) = E {log 2 (l -f SNRz^)} bits/block, (43) 

B—yO ® ^' 

i=l 


and 

ar* (SNR, 0) 


89 


^ m 

= - - C0V{log2 (1 -f SNRZ,), log2 (1 -f SNRZj ) } 


9=0 ij=l 

c 


- ^ y]E{log 2 (l + SNRZ,)} 


(44) 


where ( is dehned as 


c = 


2/3 


a(a + j5) 


(45) 
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Fig. 9. Maximum average arrival rate vs. signal-to-noise ratio SNR for 
different values of 0 and source statistics. No fading correlation, i.e., p = 0. 


Proof: See Appendix |E] 

We note that when 9 > 0, Cg depends on the source and 
channel statistics. In (l44l) . we observe the dependence of even 
the first derivative on channel correlations and source statistics 
via the covariance function and the parameter respectively. 

Next we present a high-SNR characterization of the throughput 
for Rayleigh fading. 

Theorem 6: Assume that the channel fading coefficients are 
i.i.d. in each block and fading power z = is exponentially 
distributed with unit mean (i.e., Rayleigh fading is experienced). 
Then, we have 

1 f fl^l0g2SNR + O(l) 

— r*^g(SNR, 0) = ^ PoN log 2 SNR + 0(1) if 0 < 0 < 

™ [ log2 SNR + 0(1) if 0 = 0 

(46) 

as SNR —>■ 00 . 

The proof of Theorem |6] is omitted due to its similarity to the 
proof of Theorem |3] in Appendix Similar conclusions as in 
Section III-Cll immediately apply. 

Note that the throughput expression in (l42l l suggests that for 
sufficiently high SNR levels leading to 9Ce (snr, 0) ^ a + /3, 
we have 

r*^g(SNR, 0) « PonO£;(SNR, 0). (47) 

Hence, at high SNRs, the maximum average arrival rate depends 
on the source statistics only through the ON probability. This is 
noted in the high-SNR behavior in ( |46] | as well. 

In Fig. |9] we plot r*^g vs. SNR curves for different a, /3, and 
0 values. We immediately observe that throughput diminishes 
with increasing 0 and decreasing Pqn- In Fig- [ini we analyze 
the effect of r^yg, Pqn, and a + /3 on the required SNR levels. 
For Markov fluid sources, ON state probability is not the sole 
indicator of burstiness. Having low a and /3 values also indicates 
that source is more bursty as the transition between ON and 
OFF states becomes less frequent. Hence, OFF state can be more 
persistent. When a and (3 are large, state transitions occur more 



Fig. 10. Required SNR vs. ON probability, Pqn. for a given average arrival 
rate. 0 = 0.5. No fading correlation, i.e., p = 0. 



Fig. 11. Maximum average anival rate r*yg vs. QoS exponent 6 for different 
values of a. and /3. Pqn = 0.5, p = 0, and SNR = 0 dB. 


rapidly, leading to lower required SNR levels. Again, we notice 
that the burstiness is harmful for the system. 

In Fig. [TT] we plot the maximum average arrival rate r*^,g as 
a function of 0 for different values of a and /3. Notice that by 
keeping a = (3, the ON probability Pqn is fixed at 0.5, while 
average durations of ON and OFF states vary as the values of 
a = (3 change. For example, higher a and /? values lead to shorter 
periods for ON and OFF states on average. As an outcome of this 
fact, we observe in the figure that higher throughput is achieved 
with sources having higher a + /3. 

C. Markov Modulated Poisson Sources 

In this section, we address two-state (ON/OFF) MMPP 
sources. Similarly as for the previous source models, we deter¬ 
mine the maximum average arrival rate of the MMPP source, 
which can be supported by the fading channel in the presence of 
QoS constraints, and investigate the throughput in the low-0 and 
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high-SNR regimes. The results can be immediately specialized 
to pure Poisson sources by setting /3 = 0. 

Theorem 7: For the two-state (ON/OFF) MMPP source model, 
the maximum average arrival rate is 


r*^g(SNR, e) = PoN 


9 [0Ce(snr, 9) + a + (i] 
(e®—1) [^(^^(sNR, 6) + a] 


C'_e(snr,6>). (48) 


Proof: See Appendix |F] 

It is interesting to observe that the throughput with the MMPP 
source is almost identical to that with the Markov fluid source 
model, save only for the multiplicative factor in ( |48] |. Note 
that < 1 for 9 > 0 and diminishes exponentially fast with 
increasing 9. Hence, the throughput is generally smaller with 
MMPP sources and decreases fast with 9. This can be attributed 
to the much more randomness/burstiness we experience with an 
MMPP source with respect to the previous Markov models. Note 
that the arrival rate in the ON state, rather than being a constant 
as in the previous cases, is determined by a Poisson process. 
Hence, the presence of the term is due to this Poisson 

property. Indeed, if we have a pure Poisson source, the maximum 
average arrival rate is r*^,g(SNR, 9) — (^^(SNR, 9) obtained 

by setting /3 = 0. The cost of this additional randomness is 
reflected in the following results as well. 

Theorem 8: As the statistical queueing constraints are relaxed 
by letting the QoS exponent 9 approach zero, we have 


m 

lim r*^g(SNR, 0) E {log 2 (l-f SNRZi)} bits/block, (49) 

i=l 


and 


ar* (SNR, 9) 


d9 


0=0 


= - - ^ COv{log2 (1 -f SNRZi ), log2 (1 -f SNRZj ) } 


4i=i 


^E{l0g2(l-f SNRZ,)} 




- 9 5Z®^{log2(l + SNRz0} 


2=1 


(50) 


where 


C= . 

a{a + /3) 


(51) 


Proof: See Appendix iGl 

When the system is free of QoS limitations, the maximum 
average arrival rate for the MMPP source again turns out to be 
equal to the ergodic capacity. However, the throughput has a 
steeper decline in the low-6* regime due to the third term on the 
right-hand side of JSOl l. 

Theorem 9: Assume that the channel fading coefficients are 
i.i.d. in each block and fading power z = \h\’^ is exponentially 
distributed with unit mean (i.e., Rayleigh fading is experienced). 
Then, we have 


Pq^ 


■ log2 SNR 4-0(1) if6»> 


-Cg(SNR,6») = 


log2 

^PON l0g2 SNR -f 0(1) 
Il0g2 SNR 4-0(1) 


log2 

if 0< 9< 
if 9= 0 


log2 


(52) 


as SNR —>■ oo. 



Fig. 12. Maximum average arrival rate vs. signal-to-noise ratio SNR for 
different values of 6 and different source statistics. No fading correlation, i.e., 
p = 0. 



Fig. 13. Maximum average anlval rate r*yg vs. QoS exponent 6 for different 
source statistics, p = 0. 


Since the ratio between the MMPP throughput and Markov 
fluid throughput always stays at we can immediately obtain 
the above high-SNR characterization, using the formulations in 

dSl. 

In the numerical results, we have similar conclusions as in 
the Markov fluid case. The primary difference is the reduced 
throughput for given 9, which, for instance, is readily seen when 
we compare Figs. Inland fT2l where we have throughput vs. SNR 
curves for Markov fluid and MMPP sources, respectively. In 
Fig. [El we display the maximum average arrival rate as 
a function of 9. We set a 4- /3 = 100 and SNR = 1, and vary 
a and (3 and hence the ON probability. We note that as Pqn 
decreases, the performance degrades faster with increasing 9, as 
indicated by the steeper slopes. 
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D. Comparative View of Source Models and Performance Levels 

In our analysis, we have considered discrete-time Markov, 
Markov fluid, and MMPP arrival models. All models possess 
the Markovian property in the sense that the evolution of the 
Markov chains and hence the state transitions satisfy the Markov 
condition and are described by the transition probability matrix 
in the case of discrete-time models and by the transition rate 
matrix in the case of fluid (or equivalently continuous-time) 
models. Also, state holding times are geometrically distributed in 
discrete-time models and exponentially distributed in continuous¬ 
time models, and hence exhibit the memoryless property. 

At the same time, there are distinct differences between 
different source models. For instance, transitions between states 
occur in discrete time steps in discrete-time Markov models while 
the Markov chain can spend a continuous amount of time in any 
state in Markov fluid models (i.e., the length of time spent in any 
state is a continuous random variable or more explicitly holding 
times are exponentially distributed as also noted above). MMPP 
models are further differentiated. In the discrete-time Markov and 
Markov fluid models, arrival rates are assumed to be constant in 
any given state. On the other hand, when the arrivals are modeled 
as MMPP, arrival rate is Poisson distributed in each state with a 
different intensity. Hence, MMPP sources exhibit a higher level 
of variation in this sense and can be regarded as a more bursty 
source. 

We also remark that ON/OFF discrete-time Markov and 
Markov fluid source models can be easily specialized to the 
source with a constant arrival rate by letting ON state probability 
PoN — 1- On the other hand, when Pqn = 1 in the ON/OFF 
MMPP source, we have a pure Poisson arrival source. 

Finally, we note that although there is a certain degree of 
similarity in the analysis of discrete-time Markov and Markov 
fluid models and their throughput performances (e.g., high-SNR 
characterizations are the same in Theorems [3] and |6]), the set 
of results for one model do not immediately follow from those 
for the other model as seen in the throughput formulations in 
(l34l i and (l42l) and the definitions of rj and in dJTl) and (l45T l. 
respectively. However, there is a clear distinction when MMPP 
sources are considered. As also discussed in Section ITlI-Cl higher 
level of burstiness of MMPP sources penalizes the performance, 
and lower throughput levels are achieved in general with these 
sources. 


IV. Energy Efficiency Analysis 
I n this section, we conduct a low-SNR analysis and investigate 
the energy efficiency in fading channels when data arrivals are 
random and statistical queueing constraints are imposed. We 
first identify the energy efficiency metrics. Subsequently, we 
consider different source arrival models and provide closed-form 
expressions for the energy efficiency metrics when the arrival 
rate is constant or follows a two-state Markovian model. We 
also numerically analyze specific n-state Markovian sources. 
Similarly as in the previous section, arbitrary fading correlation 
within each fading block is considered in the analysis. 

A. Energy Efficiency Metrics 

Before defining the energy efficiency metrics, we briefly de¬ 
scribe the concavity of the maximum average arrival rate as 


a function of SNR in the two-state (ON/OEE) arrival models 
(or if the arrival rates in an n-state model can be expressed as 
multiples of a certain single rate). In ll^ Lemma 1], it was 
proven that effective capacity is a concave function of SNR. 
Elwalid and Mitra m showed that the effective bandwidth of 
a source is monotonically increasing when any arrival rate Ai 
increases and is convex in the arrival rates {Ai, A 2 ,..., Aat}. 
In the ON/OEE arrival models, we have a single arrival rate 
A. Since effective bandwidth is a monotonically increasing and 
convex function of A, the inverse function of the effective 
bandwidth a*“^ exists and is a nondecreasing concave function. 
More specifically, the maximum arrival rate can be expressed 
as A*(snr, 0) = a*“^ ((^^(SNR, 6*)), which is a nondecreasing 
concave function of the effective capacity, which is concave in 
SNR. Using the composition properties of concave functions 
Eol, we realize that the maximum arrival rate is concave in 
SNR. Thus, the maximum average arrival rate r*^,g(SNR, 0) is also 
concave in SNR. 

In our analysis, following the approach in HTll . we study the 
minimum energy per bit and the wideband slope, which is defined 
as the slope of the spectral efficiency curve at zero spectral 
efficiency, as the performance metrics of energy efficiency. While 
minimum bit energy is a performance measure in the limit as 
SNR —I 0 (due to the concavity of the throughput), wideband 
slope has emerged as a tool that enables us to analyze the energy 
efficiency at low but nonzero SNR levels. In our setup, we define 
energy per bit as 


Eb _ SNR 

No ~ rJ^g(SNR,6»)/TO 


(53) 


where the normalization with m is due to our assumption that 
r*^g is in the units of bits per m symbols (or equivalently per 
block). 

The minimum energy per bit ^ . under QoS constraints can 
be obtained from 


Eb SNR 1 

— = lim -;--— = ---—. 

Wo min SNR^O r*^g(SNR, 6»)/m /•*„g(0)/m 


(54) 


At ^ , the slope So of the throughput versus Eb/No (in dB) 

curve IS defined as ED 


So = hm 

Ek. I Ek 
«0 No min 


r* (snr,6»)/to 


lOlogio - lOlogio 


Ek 

No min 


lOlogio 2 . (55) 


The wideband slope can also be found from 

C _ ^(^avg(O)/''^) , 2 


(56) 


where f*^,g(0) and f*^,g(0) are the first and second derivatives, 
respectively, of the function r*^,g(0) with respect to SNR at zero 
SNR. ^ . and iSq essentially provide a linear approximation 
of the throughput curve at low SNR levels. 


B. Energy Efficiency with Constant Arrival Rate 

In this section, we assume that the source arrival rate is fixed. 
Hence, we investigate the energy efficiency in the absence of 
source randomness and examine the impact of fading correlation 
and queueing constraints. As discussed in the previous section. 
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effective capacity, (7£;(SNR, 0), characterizes the maximum con¬ 
stant arrival rate in the presence of QoS constraints described by 
the QoS exponent 6. Hence, we in this case have r*^g(SNR, 9) = 
C£;(snr, 9). In the following result, we provide the minimum bit 
energy and wideband slope expressions under these assumptions. 

Theorem 10: Assume that the source arrival rate is constant. 
Then, the minimum energy per bit and wideband slope expres¬ 
sions as a function of the QoS exponent 9 are given, respectively, 
by 


Eh _ loge2 
TVomin E{z} 


(57) 


and 


5o = 


2(E{4)^ 


(58) 


where cov(zi, Zj) = E{ziZj} — E{zi}E{zj} is the covariance of 
Zi and Zj. 

Proof: See Appendix iHl 

Remark 2: As can be seen in ( |57] ). the minimum energy per 
bit, which is achieved in the asymptotic regime in which SNR 
vanishes, does not depend on the QoS exponent 9, hence is not 
affected by the presence of the buffer limitations. Indeed, this 
is the fundamental limit in Gaussian channels ED. Wideband 
slope 5o, on the other hand, depends on the QoS constraints 
via the QoS exponent 9. It can be easily seen that higher the 
value of 9, the stricter the QoS constraints are and the smaller 
the value of the wideband slope is, indicating the increased 
energy requirements. Furthermore, it can be readily verified that 
wideband slope decreases with increased fading correlation. Or 
conversely, variations in the channel conditions are favorable for 
improved energy efficiency. 

In Fig. [141 we plot the normalized maximum average arrival 
rate as a function of the energy per bit ^ for different 

correlation factors p when 9=1 and E{z} = E{|/ip} = cr^ = 1. 
As predicted by Theorem (TO] all curves converge to the same 
minimum energy per bit of ^ “ loge 2 = —1-59 

dB as SNR and hence r* vanish. On the other hand, wideband 
slopes are different for different values of p. As discussed above, 
as p and hence correlation diminishes from 1 to 0, slopes increase 
progressively. It is also interesting to note that in the absence of 
QoS constraints, i.e., when 0 = 0, such a distinction disappears. 
The wideband slope becomes Sq = ^ which clearly does 

not depend on the fading correlation. 


C. Energy Efficiency with Discrete-Time Markov Sources 

Starting with this subsection, we incorporate random arrivals 
into our energy efficiency analysis and determine how source 
randomness affects the performance. 

1) ON-OFF Discrete-Time Markov Sources: We assume that 
data arrival is either in the ON or OFF state in each block duration 
of m symbols. As we have previously stated, in the ON state, A 
bits arrive (i.e., the arrival rate is A bits/block) while there are 
no arrivals in the OFF state. Below, we provide our results on 
energy efficiency. 

Theorem 11: Assume that the source arrival rate is random 
and follows the described discrete-time ON-OFF Markov model. 



Fig. 14. Normalized maximum average arrival rate vs. energy per bit 

^ for different fading correlations. 9 = 1, m = 10. 


Then, the minimum energy per bit and wideband slope expres¬ 
sions as a function of the QoS exponent 9 are given, respectively, 
by 


Eh ^ loge2 
TV^Omin ]E{z} 


(59) 


and 


5o = 


2(E{z})^ 


9m 

I log^ 2 


(E{4)^ 


9 

m logg 2 


Er,, = lCOv{Zi,Zy} -f E{z 2 } 


(60) 


where rj is defined in dJTl) . 

Proof: See Appendix |l] 

Interestingly, ^ . again turns out to be a very robust quantity. 
Regardless of the "finffer constraints, channel correlations, and 
randomness of the arrivals, the minimum received energy per bit 
is = logp 2 = —1.59 dB when Ejzj = 1. On the other 

hand, the impact of random arrivals on the wideband slope is 
perspicuous in dhOl l. When compared to (fSSl) . we immediately 
notice that having random arrivals leads to the introduction of 
the term 77j^^^(E {z})^ in the denominator of (l60l) . Notice that 
when P 22 = 1 and pn = 0, we have Pqn = 1, meaning that 
we have a constant arrival rate. In this case, p = 0 and indeed 
dhOl l specializes to (l58] l. More generally, we have p > 0 for all 
Pii,P 22 G [0,1]. Therefore, random arrivals potentially decreases 
the wideband slope and increases the energy requirements. 

This is more clearly seen again in the special case in which 
pii = 1 — s and P 22 = s. Now, we have rj = ^ and the 

wideband slope is 


5o = 


1 — s 9m 
s log^ 2 


(E{z})2 


2(E{4)^ 

(61) 


As PoN = s decreases, the wideband slope decreases as well. 
Therefore, the source becoming more bursty leads to increased 
energy per bit. This is illustrated in Fig. [15] where maximum 
average arrival rate vs. energy per bit is plotted and the same 
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Fig. 15. Maximum average arrival rate energy per bit for different 

values of Pqn = ^ when p = 0.75, 9 = 1, and m = 10. 


Fig. 16. Maximum average anlval rate :^?’avg bit for different 

values of s when channel blocks are uncorrelated and 9 = 1. Number of states 
for the arrival process is n = 10. 


channel fading and correlation model as in Fig. [14] is used. In 
this figure, we assume 9 = 1 and p = 0.75. As predicted, 
the minimum bit energies are all the same. However, we have 
diminishing slopes with decreasing Pqn- Note that for a fixed 
average arrival rate, as Pqn gets smaller, source becomes more 
bursty. Data arrives less frequently but with a higher rate. This 
in turn increases energy per bit as implied by smaller wideband 
slopes. 

2) Discrete-Time Markov Sources with n States: In this model, 
we assume that there are n — 1 sources, each having its own 
ON and OFF states. In the ON state, a source sends data to the 
buffer at the rate of A bits/block. Otherwise, it is in OFF state 
in which no data is generated. In this set-up, depending on how 
many sources are active (i.e., are in ON state), data arrivals to the 
buffer can be regarded as a discrete-time Markov process with n 
states. In the state of this model, {i — 1) sources are active. 
For simplicity, we assume that the probability of each source 
being active in a given block is s, independent of the previous 
states and of the other sources. Then, the state probabilities will 
be given by 


n —1 
i — 1 


^(1 —s)" * for j = 1, 2,..., n. 


TTl 

772 • • 

• 77, 

TTl 

772 • • 

• 77, 

TTl 

772 • • 

■ 77, 


J = 


Using (l62l l we can write the average rate expression as 


I avg 



(62) 


Note that the system is essentially memoryless because each state 
is independent of the previous state. Hence, transition probability 
matrix becomes 


(63) 


= E ( ”_ - «)”■*(* - 1)^ = (" - (64) 


For this case, we do not have closed-form expressions. However, 
we can easily obtain the effective bandwidth and maximum 


Fig. 17. Maximum average arrival rate vs. energy per bit for different 

6 values when channel blocks are uncorrelated and s = 0.5. Number of states 
for anival process is n = 10. 


average arrival rate numerically. In particular, by numerically 
solving (l3^ and using (jhUi, we can determine the maximum 
average arrival rate r*^,g as a function of SNR. In Fig. [T6| we 
display the maximum average arrival rate as a function of energy 
per bit. Similarly as in the simple ON/OFF model, we observe 
that when source burstiness is decreased by increasing s, energy 
efficiency improves. 

To have a better understanding of the effect of the QoS 
constraints, average arrival rate curves in Fig. (TT] are obtained for 
different 9 values. We first notice that QoS exponent 9 does not 
have any effect on the minimum energy per bit because all curves 
merge at —1.59 dB which is again the minimum energy per bit 
for all 9 values. However, energy efficiency degrades with stricter 
QoS conditions as increasing 9 reduces the wideband slope. 

Finally, for comparison purposes, we depict the throughput as 
a function of SNR in Fig. [18] for different source characteristics 
and 9 values. The trends in the throughput vs. SNR curves for 
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Fig. 18. Maximum average aiTival rate r* vs. signal-to-noise ratio SNR when Fig. 19. Maximum average arrival rate vs. energy per bit ^ when 

channel blocks are uncorrelated. Number of states of the arrival process is n = p = 0.75, 9 = 1, and m = 10. 

10 . 


the considered n-state discrete Markov source are observed to 
be similar to those in Fig. [3] plotted for the ON-OFF discrete 
Markov source. For instance, again increased source burstiness 
(i.e., lower values of s) and stricter queueing constraints (i.e., 
higher 0 values) result in the degradation of the throughput for 
both two-state (ON/OFF) and n-state source models. 


D. Energy Efficiency with Markov Fluid Sources 

1) ON-OFF Markov Fluid Model: Now, we consider Markov 
fluid sources with two states, namely OFF state with no arrivals 
and ON state in which the andval rate is A. The generating matrix 
is defined in (l24l l. Minimum energy per bit and wideband slope 
are derived in the following result. 

Theorem 12: Assume that the source arrival is modeled by a 
two-state (ON-OFF) continuous-time Markov chain. Then, the 
minimum energy per bit and wideband slope expressions as a 
function of the QoS exponent 9 are given, respectively, by 

^ ^ l0ge2 

NoMn E{4 


and 


5o = 


2(E{4)2 


Ci^(E{4)^ + 


9 


m logg 2 




( 66 ) 


where C, = as defined in (05]). 

Proof: See Appendix U 

Similarly as before, QoS constraints and source randomness 
do not affect the minimum energy per bit. On the other hand, it 
is seen in (I66I) that the impact of source arrival characteristics on 
the wideband slope is via the state transition rates a and /3. For 
instance, larger the a value, the higher the wideband slope is. This 
is due to the fact that as a, which is the transition rate from OFF 
state to ON state, increases, the system is more likely to be in the 
ON state. Contrarily, wideband slope diminishes with increasing 


13. This is expected as well since larger (3 leads to higher OFF- 
state probabilities. The effect of a and (3 is illustrated in Fig. 
[T9| where maximum average arrival rate vs. energy per bit is 
plotted. In this figure, we set 0 = 1 and p = 0.75. As predicted, 
the same minimum energy per bit is achieved for different values 
of a and /3, while wideband slope increases with increasing a 
or decreasing f3. 

2) n-State Markov Fluid Birth-Death Process: In this sub¬ 
section, we consider a birth-death process for the Markov fluid 
source. We assume that there are n states and the arrival rate in 
the state is (r — 1)A for i = 1,..., n. The generating matrix 
for the birth-death process is in the form of 


—a a 0 • • ■ 

f3 —(a-l-/3) a 0 



0 ••• 0/3 

0 . 0 


O' 

0 


0 

— (a-|-/3) a 
-/3 /3. 


(67) 


Hence, the transition rate from state i to state z -|-1 is a whereas 
the transition rate from state i to state z — 1 is /3. The effective 
bandwidth of this source, which does not have a simple closed- 
form expression, can be found from (1221) . In order to conduct 
an energy efficiency analysis, average arrival rate needs to be 
identified as well. Using (|2^ and dhTl) . we can easily determine 
that the stationary distribution as 


Til = 


(?) 

i-1 

1 

- 1 

(?)■ 

1 

- 1 

(?) 

n 

1 


for z = 1, 2,..., rz 


( 68 ) 


when j 1, a 0 and /3 7 ^ 0. If any of these inequalities 
is not satisfied, state probabilities can be obtained by limiting 
functions. 

Now, under the assumptions that a /3 and the arrival rate in 
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Fig. 20. Maximum average arrival rate :^’"avg energy per bit when 
channel blocks are uncon'elated. 6 = 1, f3 = 100 and the number of states of 
the arrival process is n = 10. 



Fig. 21. Maximum average arrival rate r^vg vs. energy per bit when channel 
blocks are uncorrelated and a = /3 = 50. Number of states of the arrival model 
is 10. 


State OS = (i — I)A(snr), the average arrival rate is given by and 




(69) 


where ^ 

Remark 3: Note that ( |69] | specializes to (l27l l if n = 2. When 
^ ^ oo, the probability of the state approaches 1 , and the 
arrival rate will be (n — I)A(snr) at steady state. On the other 
hand, for ^ = 0 , the state of the source is stuck at the first state 
in which the arrival rate is zero. 

Numerically, we can obtain the effective bandwidth of the n- 
state birth-death Markov fluid process using (l22T l. Subsequently, 
solving (l3?t and incorporating ( l69b , we can determine the maxi¬ 
mum average arrival rate r*^,g(SNR), which we further employ for 
characterizing the energy efficiency. The results of this numerical 
analysis are displayed in the following figures. In Fig. |20] we 
demonstrate the effect of a on the energy efficiency. In particular, 
when P is kept fixed, increasing a improves the energy efficiency 
as in the two-state case. We illustrate the effect of QoS constraints 
in Fig. |2T] Similar conclusions as before readily apply. QoS 
exponent 9 does not alter the minimum bit energy, which is —1.59 
dB again, but the wideband slope is reduced with increasing 9. 


E. Energy Efficiency with Markov-Modulated Poisson Process 

1) ON-OFF Markov-Modulated Poisson Process: Again, we 
initially address the two-state model in which there are no arrivals 
in the OFF state and the intensity of the Poisson arrival process 
is A in the ON state. The generating matrix G is the same as in 

(EUl. 

Theorem 13: Assume that the source arrival is modeled by a 
two-state (ON/OFF) Markov modulated Poisson process. Then, 
the minimum energy per bit and wideband slope expressions as 
a function of the QoS exponent 9 are given, respectively, by 


^ ^ - 1 loge 2 

V^Omin S E{z} 


(70) 


5o = 


Ci^(E{z})^ 


(71) 


where ( = as defined in ( |45] |. 

Proof: See Appendix IkI 

Remark 4: It is interesting to observe that, unlike the previous 
arrival models, minimum energy per bit in the case of MMPP 
source depends on the QoS exponent 9. More specifically, 
minimum energy per bit increases with (e® — l)/9 which is an 
increasing monotonic function of 9 and always greater than one 
for 0 > 0. On the other hand, as 9^0, (e® — l)/0 —^ 1. There¬ 
fore, ^ with equality only if no QoS constraints 

are imposed (i.e., when 9 = 0). Furthermore, in addition to its 
significant impact on the minimum energy per bit, increasing 9 
leads to much quicker reduction in the wideband slope due to 
the presence of the term in GD- Hence, overall, energy 

costs grow very fast as 9 increases. This is again because of the 
additional randomness arising from Poisson arrivals in the ON 
state. 

Remark 5: From dZB, we note that the effect of the state 
transition rates a and /3 on the energy efficiency is the same as 
in the Markov fluid source model. Increasing a or decreasing 
P improves the energy efficiency of the system because the 
burstiness of the data arrivals is reduced and the buffer overflows 
can be avoided at lower energy costs. 

We plot the maximum average arrival rate vs. energy per bit 
in Fig. |22] We set Ejz} = 1 and 9 = 1 for which the minimum 
energy per bit is 0.76 dB. The increase in bit energy with respect 
to —1.59 dB is due to 101 og 2 Q((e® — l)/9) = 101ogig(e — 1) 
for 9 = 1. From the figure, we can again infer that adjusting a 
or P to increase the ON state probability makes the system more 
energy efficient due to the increase in the wideband slope. 

2) n-State Markov-Modulated Poisson Process: Finally, we 
consider an n-state MMPP process and assume that the intensity 
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Fig. 22. Maximum average arrival rate energy per bit when 

p = 0.75, 0=1, and m = 10. 


of the Poisson arrivals in the state is (i — 1)A. For the Markov 
transitions between states, we consider the birth-death process 
and adopt the transition rate matrix G from dhTl i. We solve for the 
maximum intensity A*(SNR, 0) by incorporating (l28T l into (l3Al l. 
Then, using the expression in (|69] |, we obtain r*^g(SNR, 0). 

In Figs. 122 and |24] we depict the maximum average arrival 
rate as a function of the energy per bit with uncorrelated channel 
coefficients being assumed in each block. In Fig. |23] we set n = 
10 , 0 = 1 and (3 = 100, and demonstrate how a influences 
the energy efficiency of system. The observation has similarities 
with other Markovian sources regarding the source burstiness. 
Interestingly, the minimum energy per bit is again 0.76 dB as 
in the two-state case, leading to the conclusion that the number 
of states does not alter ^ in this case. The degradation in 

-^Omm ® I_I 

energy efficiency due to increased 9 is shown in Fig. [2^ As 
described in the two-state case, higher values of 9 (i.e., stricter 
QoS constraints) result in higher ^ . and smaller wideband 
slope. Therefore, even for relatively small increases in 9, we 
can have large gaps between curves, indicating significantly high 
energy costs. 

V. Conclusion 

In this paper, we have studied the throughput and energy 
efficiency in wireless fading channels when data arrivals to the 
buffer are random and constraints on buffer overflow probabilities 
are imposed. We have considered discrete-time Markov, Markov 
fluid, and Markov-modulated Poisson (MMPP) sources, and 
formulated the maximum average arrival rates of these sources 
in closed-form in the special case in which the source has only 
ON and OFF states. We have shown that as the QoS exponent 
9 vanishes (i.e., no QoS constraints are imposed), throughput 
converges to the ergodic capacity and the performance becomes 
independent of the source characteristics and channel fading 
correlations. On the other hand, we have demonstrated for 9 > Q 
that the statistics of the source arrivals and channel conditions 
have significant impact on the throughput. More specifically, it 
is seen that while fast channel variations (e.g., having fading less 



Fig. 23. Maximum average arrival rate vs. energy per bit when 

channel blocks are uncorrelated. 6 = 1, fi = 100 and the number of *^arkov 
states for the arrival process is n = 10. 



Fig. 24. Maximum average arrival rate vs. energy per bit when 

channel blocks are uncorrelated and a = fi = 50. Number of Markov states for 
the arrival process is n = 10. 


correlated) can be beneficial, source randomness and burstiness 
degrade the throughput. In particular, performance loss can be 
severe with MMPP sources especially if 9 is large. Similar con¬ 
clusions are drawn in the energy efficiency analysis as well. We 
have shown that the minimum energy per bit is only a function 
of the marginal distribution of the fading for constant, discrete¬ 
time Markov and Markov fluid arrivals. However, for MMPP 
arrivals, this asymptotic measure of energy efficiency increases 
as the QoS exponent 9 increases and hence as stricter queueing 
constraints are imposed. For practical scenarios, wideband slope 
is a better indicator of the performance and is always seen to 
depend on the source statistics, channel correlations, and QoS 
exponent 9. In particular, we have noted that wideband slope 
diminishes with increasing source burstiness and stricter QoS 
constraints, indicating increased energy costs in such cases. While 
this paper primarily concentrates on two-state Markovian sources, 
numerical results are provided in the energy efficiency analysis 
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for several simple n-state models with which we have similar 
observations as above. 


Appendix 

A. Proof of Theorem Q} 

Using the effective bandwidth formulation in (fT^ . we can 
express (l3^ in the following equivalent form; 



pil+P22e''® + -y/(pI7+P22e^^®U--4(pii+P22--l)e^ 
2 


Ce. 


(72) 


Then, we rewrite the above equality as 

Pit +P22e^® + \J{Pii + P22e^®)2 - 4(pii + P22 - l)e^® = 

(73) 


from which, after moving the hrst two terms on the left-hand 
side to the right-hand side and taking the square of both sides, 
we obtain 


(Pii + P22e^®)^-4(pii-f p22-l)e^® = (2e®'^® - Pii-P22e^^^. 

(74) 


Now, by simply exchanging the second term on the left-hand side 
with the term on the right-hand side, we have 

(Pii + P22e^®)^- (2e®'^® -Pii-P22e^^^= 4(pii-f p22-l)e^®, 

(75) 

(2pn+ 2 p 22 e^®- 2e®^®) 2e®^® = 4(pn + p22 - l)e^^ 

(76) 

After further rearrangements, we have 

(Pii +P 22 - 1 -P 22 e®^®)e^® = (77) 


Solving the equation for A, we get 


A*(snr, 0) = ilog. 


g2ecE(SNR,e) _p^^gecE(SNR,e) 

1 - Pit - P 22 + ^ 226 ®*^®*^™^’®) 


(78) 


which provides the maximum arrival rate in the ON state. We 
can now express the maximum arrival rate as r* (SNR, 6) = 
-PonA*(snr, 0) and obtain the expression in (l34l) . □ 


B. Proof of Theorem \2j; 

Let us dehne 

i;{e) = e-«f^®(SNR ,0 . 

(79) 

The following properties of t/j can be verihed easily: 

i/’(0) =1, (80) 

'f{9) =E | — log 2 (l + SNRZi)e~^^‘=i l°g 2 (l+SNRzi) I ^ 


i=l 


(81) 


IA( 6 >) =E I ^log2(l -f SNRZ,) 


(82) 


and 


lp{0) = - E l0g2(l -f SNRZj) !> , 

■ip{0) =E<| ( ^log 2 (l + SNRZi) ] 




(83) 

(84) 


where tp and ip denote the hrst and second derivatives of ip with 
respect to 9, respectively. Additionally, we dehne (SNR, 9) as 

Cg(SNR,6l) = (85) 

Therefore, by applying L’Hopital’s rule and letting 6 * —0, 
maximum average arrival rate and its slope can be easily found 
as 


limr* (SNR, 0)= /i(0), 

gCg(SNR,0) ^ /i(0) 

39 2 


( 86 ) 

(87) 


Now, replacing e ®C'e(SNR. 6 ») jn expression of 

r*^g(SNR,0) in (1^ , we can express fi{9) as 


/l((9)=PoN[logg(l-pill/’(^'))-loge((l-Pll-P22)V’^(6')+P22V'(6’))] • 

( 88 ) 

Therefore, we derive fi{9) when 0 —0 as 


/l(0) =P0N 


-P11'0(O) 

l-pil'0(O) 


(l-pil-p22)l/’(0) 

(1 -pil -P22)V'(0) +P 22 


= PoN 


-Pit 


.1-Pll 

^ 771 

= Ei ^log 2 (l 


1 -P 11 -P 22 

1-pll 
■f SNRzU > . 


- 1 


^( 0 ) 


. i=l 


^( 0 ) 

V'(o) 

(89) 

(90) 

(91) 


Note that (l89l l follows by taking the hrst derivative of the 
expression in dHSl) with respect to 9, and (|90] | is obtained using 
the property that ip{0) = 1. Finally, ( 1911 1 and hence the result in 
drSl l immediately follow from ( l20l l, ( l83l l and dS^ . 

Next, we determine the slope of the throughput in (l36l) as the 
QoS exponent 9 approaches zero. For this, we only need to derive 
the second derivative expression /i( 0 ), which is done at the top 
of the next page. (|92] |. (l93l l and (l94l l follow from straightforward 
algebraic steps. Inserting (ISTl l and (l84l l into (l94l l. we obtain ( |95]) . 
Finally, the result in (l36l l follows by combining (l95l l and dSTl i. □ 


C. Proof of Theorem |2} 

In the analysis of the high-SNR slope of the effective capacity, 
it has been shown in 1281 that 


- i logg E |e“® log 2 (l+SNR 2 ;) I 


gj^l 0 g 2 SNR + C9(l) if 0 > 
log2 SNR-1-0(1) 


log 2 e 

if 0 < 6 » < 1 


(96) 


log2 e 
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/.( 0 ) = F„.^ 


= -PoN 


-P110(O) 

(l-Pll-P 22 ) 0 (O) 

0 ( 0)1 

1-P110(O) 

(1 -Pll -^ 22 ) 0 ( 0 ) +P 22 

0(O)J 


-Pll I-P 11 -P 22 


I-Pll 


I-Pll 


-1 


^( 0 ) 


9=0 

Pll , (1-Pll-P22)^ ' 

(l-pn)2 ^ (l-pii)2 ^ . 


= -^(0) + (1 - ?y) (^(0 


= -E. 


^l0g2(l + SNRZi) ] I +(1-77) 


IE< ^l0g2(l +SNRZi) 



(92) 

(93) 

(94) 

(95) 


where z is exponentially distributed with unit mean. If we assume 
that fading in each block is i.i.d., then the effective capacity 
expression in (fOl l becomes 

CE{SNR,e) = -ilog^E{e-®^r=ii°g(i+SNR^O| (97) 

= -Jloge (98) 

=log^E{e-®^°«(^+™'^^)}. (99) 

Furthermore, the maximum average arrival rate in (l34l i can be 

expressed as 

r:,g(SNR, 0) 

Po^ ( e2ecE(SNR.e)(i_p^^g-ec^(SNR.e)^ \ 

- — loge(^geCE(SNR,e) ((1 _ _ P22) e-«CB(SNR.e) + P22) ) 

( 100 ) 

PoN , f (1 _p^^g-.C.(SNR,«)) 

= ^ (^log, + log, (1 - 

- log, ((1 - pn - P 22 ) + P 22 ) ) 

( 102 ) 

= ^ log, + 0(1) (103) 

0 

= PonC'b(snr,0) + O(1) (104) 

where ( llOll i and ( 11021 ) follow from straightforward algebraic 
operations and ( 11031 ) is due to the fact that Ce (snr, 0) increases 
without bound as SNR increases and hence the term e-®C'E(SNR,6i) 
vanishes asymptotically in the formulations. 

Finally, combining (l96l ), (l99l ), and (11041 ), we immediately 
obtain the desired result in (|3^ for the cases in which 0 > 0. 
When 0 = 0, the result follows from ([35]) in Theorem |2| □ 

D. Proof of Theorem [4} 

Using (1251 ). we can rewrite ([33]) as 

{0\ -{a+ 13)- 20 cEf = (0A - (a + (3))^ + Aa0\ (105) 
which can further be simplified to 

-20Ce{ 20\ - 2{a + (3)- 20 Ce) = 4a0A. (106) 


Next, solving for A, we obtain 

A* (SNR, 0) = I ^ ^ gg(SNR, 0). (107) 

0C_e(snr, 0) + a 

Finally, using the expression in (1271 ). we derive the maximum 
average arrival rate given in (|42] ). □ 


E. Proof of Theorem [5] 

Similar as in the Proof of Theorem[2]in AnnendixlBl we define 
r*^g(SNR, 0 ) = ^ with 

/ 2 ( 0 ) = log, m- (108) 

a - log, 'f{0) 

Now, we have (1861 ) and dSTl ) hold with fi replaced with / 2 . The 
remainder of the proof requires only the determination of the 
first and second derivatives of / 2 ( 0 ) at 0 = 0. The first derivative 
72 ( 0 ) is given at the top of the next page in (I1091 )- (ll 121) . Note that 
(11091 ) and (II101 ) follow from straightforward algebraic steps, and 
(111 11 ) is obtained by noting the property that '0(0) = 1. Finally, 
(II121 ) and hence the result in (|4^ immediately follow from ([26|, 
([83]) and ([ 86 ]). 

Next, we obtain the slope expression in (l44l ) in the limit as the 
QoS exponent 0 approaches zero. For this, we characterize the 
second derivative expression 02 ( 0 ) on the next page in ( II131 )- 
dl 161 ) . dl 131 ). dl 141 ) are readily obtained and dl 151 ) is determined 
by noting that '0(0) = 1. We incorporate (l83l ) and (l84l ) into dl 151 ) 
to obtain dl 161 ). The result in (l44l ) follows by combining dl 161 ) 
and (l87l ) (with /i replaced with 02 ). □ 


F. Proof of Theorem [7] 

We find the maximum average arrival rate r*.^,g(SNR, 0) by 
incorporating ( [29] ) into ([33]) and expressing (l33l ) as 

{{e^ -l)\-{a+ p)-20cEf 

= ((e®-l)A-(a + /3))" + 4a(e®-l)A. (117) 


Similarly as in the proof of Theorem|4] we can simplify the above 
equality and solve for the maximum Poisson arrival intensity in 
the ON state to obtain 


A* (snr, 0) 


^ [0C'_e(snr, 0) + a + /3] 
(e®— 1) [0C'£;(snr, 0) + a] 


C'e(snr,0). 


(118) 


With this characterization, the maximum average arrival rate is 
readily obtained from (l30l ). □ 
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/ 2 ( 0 ) = lim -PoN< 


= lim -PoN 
d->-0 


.|g) (a + /?_log^^(0))(-|gi) 


a-log^T/j(0) 

a/3 


(a - log^ipiff))^] - 1 ^( 0 ) 


(a - logg 'i/;(0))^ 

iP(0) 


loge ^(0) 


g + /3 - loge i/^(0) i){0) 
a-logeV'(fi') V'W 


= -Po^Mo) 


= Ei^log2(l + 

li=l 


SNRZi 


(109) 

( 110 ) 
( 111 ) 
( 112 ) 


MO) = lim-PoNjj 


= lim -PoN- 


1 + 


aP 


Vi(6») 


{a - \og^il){0)Y\ 'il){0) 

("Im) ('ip{0)\ , 


e-s-o dO I (a — logg'(/'(^))^ \'0(^) 


1 + 


af3 


{a - log^ip{0))‘^ 


= -m + (1 - 


2/3 


Q;(a + l3) 


) (' 0(0 


= -E < ( ^l0g2(l +SNRZ,) ] >+(l- 


2/3 


a(Q; + /3) 


0 W _ ( 

Vi( 6 i) I V'( 6 ') 


) E ^log2(l + 


SNRZ,- 


(113) 

(114) 

(115) 

(116) 


G. Proof of Theorem |S} 

Employing f2{0) defined in (1108b . we can express the maxi¬ 
mum average arrival rate as 


Cg(SNR, 0) = 


hm_ 

e® — 1 


(119) 


Then, the throughput in the limit as 0 approaches zero is given 
by 

lim Cg (SNR, 0) = liin = /2(0). (120) 

Inserting the result from (II 12b into (1120b . we obtain (l49l l. Next, 
we determine the slope of the throughput when 0 approaches 
zero; 


a-* (SNR, 0) 


d0 


= i.m 4 a 

e^o — 1 


6 ^ 2 ( 0 ) 


9=0 c - i. (e® - 1 ) 

(e^-l)/2(g)-eV2(g) 
{eO - 1 )^ 

(e^-1)/ 2 (g) - 6 ^ 2 ( 0 ) 

2 (e® - 1 ) e® 


= lim 
6>^0 


= lim 
6>^0 


_ i lini / 2 (g) 


2 2 e® - 1 

/ 2 ( 0 ) / 2 ( 0 ) 


( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 


(1121b follows by taking the derivative of the expression in (II 19b 
with respect to 0. (1122b is obtained by simplifying (1121b . We 
apply L’Hopital’s rule on (1122b to get (1123b and further simplify 
it in (1124b . Finally, we obtain (1125b . which we used to derive 
(l50b by inserting (II 12b and (II 16b into (1125b . 


H. Proof of Theorem 1701 

When the arrival rate is hxed, the following equality holds: 


r:^g(SNR, g) = (3'£;(SNR, g). (126) 


Therefore, in formulas (l5^ . (l56T l. we can use C£:(0) and C£;(0) 
instead of r*^,g(0)/TO and rl^g{Q)/rn respectively, where we have 
defined C£;(SNR, 0) = Ce {SNR, 0)/m as the normalized effective 
capacity. Minimum energy per bit and wideband slope becomes 


and 


^ ^ 1 
3V0min C£(0)’ 


5o = - 


2 (Cr,( 0 ))^ 

Cb(0) 


loge 2. 


(127) 

(128) 


Thus, we only need to obtain the hrst and second derivatives 
of Cb (snr, g) with respect to SNR at SNR = 0 to determine the 
minimum energy per bit and wideband slope. We first express 
the effective capacity given in (fTSb as 


Cb(SNR) = -_loggg(SNR) 
dm 


where we have defined 




(129) 

(130) 


Now, the hrst and second derivatives of C£;(snr) with respect to 
SNR are easily seen to be given by 


Ce(snr) = - 


1 g(SNR) 


and 


Ce(snr) = - 


0m g(SNR) ’ 
g(SNR)g(SNR) - [g(SNR)]2 


0m 


[g(SNR)]2 


(131) 

(132) 
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where g and g denote the first and second derivatives of the Next, we let SNR 0. Since the arrival rate A —^ 0 when SNR —?> 
function g with respect to SNR and can be expressed as 0, the equality in (I1381 l becomes 


g(SNR)= - 

and 

g(SNR) 


-E • ' 


loge 2 1^1+ SNRZj 


„-7;n472 E” 1 ioge(i+SNRzi 


(133) 


-A E; 


0 


E 




logg2 U (l+SNRZi)2 logg2 (l+SNRZi)(l+SNRZj) 

\\ t —i 't,J — -L / 

9 


X e 


;2 E“iloge(l+SNRzi) 


(134) 


Then, at SNR = 0, we have 

e{4 


C£;(0) = 


mloge2 loge2 


(135) 


and 

Ci5(0) 


m(logg 2)2 

^ cov {z„ Zj} + m log^E {z^} 


TO(l 0 gg 2)2 


(136) 


where we have used the facts that = rnMj{z} and 

IE {-z?} = rnE{zA due to our assumption that the fading 
coefficients and therefore {zi}’s are identically distributed. 

Plugging the expressions in (1135b and (1136b into those in (1127b 
and ( 1128b . we readily obtain the minimum energy per bit and 
wideband slope expressions in ( |57] | and ( l58l l. ■ 


I. Proof of Theorem 1771 

To show the result, we need to obtain the first and second 
derivatives of r^(SNR). We first express the maximum average 
arrival rate in (l3^ as 


avg 


(SNR, 0) = ^ 


logg(l - piig(SNR)) - logg(g(SNR)) 


- logg ((1 - pii - P22)g(SNR) + P 22 ) 

(137) 


where we have used the definition that with 

g(SNR) defined in (1130b . Taking the first derivative with respect 
to SNR, we obtain 


Cg(04) 



Pit _ ^ _ 1 -Pii -P 22 

l-pii 1-pii 

(139) 


m 

0 


logg 2 


Ee{e} 


mE{z} 

logg 2 


(140) 


where Pqn = 2 -jni-p 22 ' ’■I'® result in (1140b into (l54l) . 

we immediately obtain (l5^ . 

In order to find the wideband slope, we first determine the 
second derivative of the maximum average arrival rate with 
respect to SNR and then evaluate it at SNR = 0 as follows: 


KA0,s) = 


m 


P( 


ON 


Pit 

1 -Pll 


- 1 - 


[g(0)]^ 


p, 


ON 


Pit 


(1 -Pll -P 22 ) 

1 -pll 

(1 -Pll -P 22 Y 




2 


1 + 


0 


( 1 -Pll )2 

(141) 


(1141b follows from the fact that p(0) = 1, and p is defined in 
(IJTI i. Finally, inserting (1140b and (1141b into (l56l) . the wideband 
slope expression in (l60l) is readily obtained. ■ 


J. Proof of Theorem 1721 

We differentiate the maximum average arrival rate expression 
in (l42b with respect to SNR and obtain 


f* (SNR, 0) = PoN 


26»(7£;(snr)(7£;(snr) + (o + ,8)C'_e(snr) 


6»(7£;(snr) + a 
[6>(7|(snr) + (a + /3 )C'_e(snr)] 6>(7e(snr) 1 


(6»(7£;(snr) + a) 


r 


(142) 


As SNR — 7 > 0, we can easily derive 

4g(0, 0) = PoN^^C'i5(0) = Cs(0) = (143) 

® a logg 2 

where we use the facts that C£:( 0 ) = 0 and Pqn = 

the result in ( 1143b into ( l54b . we immediately obtain ( l65l l. 

In order to determine the wideband slope, we additionally take 
the second derivative of the maximum average arrival rate with 
respect to SNR and evaluate it at SNR = 0 as 

4v,(0, 0) = Ce( 0) - (cE(O)y. (144) 

Now, inserting the results in ( 1143b and (1144b into ( |56] ) and using 
the formulations in (1135b and (1136b . we obtain (l 66 b . ■ 


4g(SNR,6») 


Pqn -piig(SNR) 

0 l-piig(SNR) 


g(SNR) 

g(SNR) 


(1 -Pll -p22)g(SNR) 

(1 - Pll - P22)g(SNR) + P 22 


K. Proof of Theorem 1751 

The proof is rather straightforward after realizing that 
4g(SNR, 0) of the MMPP source given in (1481) is equal to the 
(138) maximum average arrival rate of the Markov fluid source in 
(l42b scaled with Therefore, making use of the results in 
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(fT43l) and dHH, we can immediately express the first and second 
derivatives of r*„g(SNR, 0) at SNR = 0 as 

0mE{z} 


K.M0) = 




(e®— 1) (e®— 1) logg 2 ’ 

Ce{ 0) - {cE{0)y 


(145) 


[24] 

[25] 

[26] 


a{a + (3) 

Then, the expressions in dTOb and ( fTTl ) are obtained by plugging 
(11451) and (11461) into (l54l) and (1561) . ■ 


(146) [27] 
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